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ABSTRACT 

Let  A  be  an  infinite  sign  regular  (SR)  matrix  which  can  be  viewed  as  a 
bounded  linear  operator  from  to  itself.  It  is  proved  here  that  if  the 

range  of  A  contains  the  sequence  (••*,  1,-1,1,-1,**»)»  then  A  is  onto. 
If  A~^  exists,  then  DA-1D  is  also  SR  ,  where  D  is  the  diagonal  matrix 
with  diagonal  entries  alternately  1  and  -1.  In  case  A  is  totally 
positive  (TP),  then  DA-1D  is  also  TP  under  additional  assumptions  on  A. 
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SIGNIFICANCE  AND  EXPLANATION 


Spline  approximation  is  often  most  effective  when  the  breakpoint  (knot) 
sequence  can  be  chosen  suitably  non-uniform.  At  the  same  time,  the  standard 
approximation  schemes  (such  as  least  squares  approximation,  or  interpolation 
at  suitable  interpolation  points  by  splines)  are  so  far  only  known  to  be 
usable  and  bounded  as  long  as  the  breakpoint  sequence  is  almost  uniform.  The 
problem  of  showing  existence  and  uniqueness  of  bounded  spline  approximants  to 
bounded  data  boils  down  to  showing  invertibility  of  a  certain  infinite 
matrix  a  .  The  distinguished  feature  of  this  matrix  is  its  total  positivity, 
i.e.,  all  minors  of  A  are  nonnegative.  In  this  paper  we  show  that  if  the 
range  of  an  infinite  totally  positive  (more  generally  sign  regular)  matrix 
A  contains  the  particular  sequence  ( • • • ,-1 , 1 ,-1 , 1 ,-1 , ♦ • • )  ,  then  every 

bounded  sequence  is  contained  in  the  range  of  A.  In  spline  terms,  this 
result  says,  for  example,  that  any  bounded  data  sequence  can  be  interpolated 
with  a  bounded  spline  (with  a  given  knot  sequence,  at  a  given  interpolation 
point  sequence)  provided  that  the  periodic  data  {+1,-1}  can  be  interpolated 
by  a  bounded  spline  from  that  class.  Our  arguments  show  that  such  an 
interpolating  spline  can  be  'constructed'  as  a  limit  point  of  splines  which 
satisfy  finitely  many  of  the  given  interpolation  conditions  provided  that  the 
trivial  data  can  be  interpolated  only  by  the  zero  spline. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


Inverses  of  infinite  sign  regular  matrices 


by  C.  de  Boor,  S.  Friedland  and  A.  Pinkus 

0.  Introduction*  If  the  problem  of  spline  interpolation  is  expressed  in  terms  of 

B-splines,  then  the  question  of  existence  of  a  bounded  spline  interpolant  to  bounded  data 

is  seen  to  be  equivalent  to  the  question  of  whether  a  certain  bounded  band  matrix  has  all 

bounded  sequences  in  its  range.  In  [4] ,  C.  A.  Micchelli  conjectured  that  there  exists  a 

unique  bounded  spline  interpolant  (of  a  given  order  and  a  given  knot  sequence)  to  any  data 

sequence  (T  ,y  ).  in  the  plane,  with  (t  )  strictly  increasing  and  (y.)  bounded, 
i  i  x€Z  1  1 

provided  only  that  it  is  possible  to  interpolate  the  particular  data  sequence 

(t  ,  ("I)1),,.,,  by  such  a  spline.  There  is  apparently  nothing  special  about  the  particular 
1  lcZ 

spline  problem  other  than  that  it  leads  to  a  banded  totally  positive  matrix.  Therefore  one 
of  us  quoted  this  conjecture  in  [2)p.319,  Problem  4)  as 

"A  biinfinite  banded  totally  positive  matrix  A  is  boundedly  invertible  if  and 
only  if  the  linear  system  Ax  »  ((-J)1)  has  a  bounded  solution." 

In  hindsight,  it  is  easy  to  see  that  this  conjecture  is  faulty  even  in  the  original  context 
of  spline  interpolation.  For  example,  interpolation  by  bounded  broken  lines  with  breakpoint 
sequence  Z  at  the  sequence  T  =  z\{0}  is  possible  to  any  bounded  ordinate  sequence  y  , 
but  not  uniquely  so  since  the  value  of  the  interpolant  at  0  is  freely  choosable.  In  ma¬ 
trix  terms,  this  corresponds  to  the  matrix  obtained  from  the  biinfinite  identity  matrix  by 
dropping  one  row.  But,  with  the  condition  changed  to  "...  has  a  unique  bounded  solution", 
the  conjecture  was  proved  in  f  1 J . 

The  argument  in  (11  establishes  that,  under  the  given  condition,  A  "has  a  main  dia¬ 
gonal",  i.e.,  some  diagonal  of  A  has  the  property  that  all  finite  sections  of  A  having 
a  portion  of  this  diagonal  as  their  main  diagonal  are  invertible,  with  their  inverse  bound- 


*)  Actually,  one  of  the  editors,  enlightened  by  [1],  changed  it  to  "a  unique 
bounded  solution"  as  the  book  went  to  press. 
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ed  uniformly.  A~^  is  obtained  as  the  pointwise  limit  of  these  inverses.  Thus,  the  argument 
establishes  more  than  Hicchelli's  conjecture.  In  reaction  to  a  presentation  of  these  argu¬ 
ments,  one  of  us  sugggested  that  there  might  be  simpler  ways  to  establish  the  conjecture 
directly.  In  particular,  it  should  be  possible,  because  of  the  checkerboard  nature  of  in¬ 
verses  of  totally  positive  matrices,  to  establish  that  A  is  onto  under  the  original 
condition,  using  minimal  solutions  of  finite  sections  of  the  given  linear  system  Ay  -  v  . 

The  present  paper  carries  out  this  program  in  Section  1.  As  it  turns  out,  it  is  pos¬ 
sible  (i)  to  drop  any  kind  of  structure  assumption  on  A  such  as  handedness,  and,  less 
surprising,  (ii)  the  assumption  of  total  positivity  can  be  relaxed  to  sign  regularity. 

Having  settled  this  matter,  it  then  became  of  interest  to  see  how  much  more  informat¬ 
ion  about  the  inverse  of  a  totally  positive  matrix  could  be  obtained  by  this  approach. 
Specifically,  assuming  A-1  to  exist,  and  with  D  the  diagonal  matrix  having  alternately 
1  and  -1  on  its  diagonal,  could  (i)  the  sign  regularity  of  DA-1D  be  established, 

(ii)  DA-1D  or  its  negative  be  shown  to  be  totally  positive  if  A  is,  (iii)  A-1  be 
approached  by  inverses  of  finite  sections  of  A  ? 

As  to  the  third  question,  we  show,  as  a  simple  corollary  to  the  development  in  Section 
1,  that  A-1  can  indeed  be  approached  pointwise  by  inverses  of  certain  submatrices  of  A 
(involving  consecutive  columns  of  A  but  not  necessarily  consecutive  rows),  provided  the 
columns  of  A  are  already  in  Cg  and  not  just  bounded.  We  believe  this  assumption  to  be 
unnecessary  in  case  A  is  totally  positive,  in  the  sense  that  we  believe  the  columns  of  a 
totally  positive  invertible  matrix  to  be  already  in  Cg  .  But  we  have  not  been  able  to 

prove  this.  In  any  case,  while  this  result  is  far  from  establishing  that  such  A  has  a 
main  diagonal,  it  does  allow  the  conclusion  that  DA_1D  or  its  negative  is  totally 
positive  in  case  A  is. 

As  to  the  first  two  questions,  we  show  in  Section  2  by  a  completely  different  line  of 
reasoning  that  DA_1D  must  again  be  sign  regular.  From  this,  a  surprisingly  simple  argu¬ 
ment  proves  the  total  positivity  of  DA_1D  in  case  A  is  totally  positive  and  infinite 
but  not  biinfinite. 

we  will  use  the  following  notations  and  conventions. 


2 


We  use  lower  case  letters  to  denote  elements  of  R1  ,  i.e.,  real  functions  (or,  se¬ 
quences)  on  some  integer  interval  I  ,  with  v(i)  the  i-th  entry,  or  value  at  i  ,  of  the 
sequence  v  ■  By  S~(v)  we  mean  the  number  of  strong  sign  changes  in  the  sequence  v  , 


(v)  sup  {  r  s  there  exist  j1  <  ...  <  jf+1  s.t.  v(j  )v(j8+1>  <  0  }  , 


S  (v)  :«  sup  {  S  (w)  :  w( i)  -  v(i)  whenever  v(i)  *  0  } 

denotes  the  weak  sign  changes  of  v  If  J  is  a  subset  of  I  ,  then  Vj  denotes  the 

restriction  of  v  to  J  while  v^j  is  shorthand  for  the  restriction  of  v  to  IS J  , 

i.e.,  to  the  complement  of  J  in  I  .  If  J  cons' sts  of  just  one  point,  J  «  {j}  say, 

then  we  write  \j  instead  of  \{j}  .  Also, 


while,  if  also  u  e  R  ,  then 


I v | ( i )  :«  (v(i) |,  all  i  , 


u  v  :=  Zigl  uli)v(i) 


Correspondingly,  when  also  J  is  an  integer  interval,  then  A  denotes  the  transpose 
Jx  J 

of  the  matrix  A  e  R  and  AR  L  denotes  the  restriction  of  A  to  the  subset  K*L  of 
I*J  .  Such  a  matrix  A  is  sign  regular  («:  SR)  provided  that  for  each  k*1,2,3,...  all 
minors  of  A  of  order  k  have  the  same  sign.  If  this  sign  is  positive  for  all  k  ,  then 
A  is  totally  positive  (“:  TP).  We  denote  the  minor  of  A  obtained  from  rows  p  <  .,.<  q 


and  columns  r  <  ...<  a  by 


a(p'  q)  . 
'r,  ...,  s' 


1.  Existence  of  a  bounded  right  inverse  in  same  absolute  norm. 

Let  J  be  a  finite,  infinite  or  biinfinite  integer  interval  and  let  S  C  If1  be  a 


norraed  linear  space  of  real  functions  on  J  ,  i.e.,  a  space  of  sequences.  We  assume  that 
the  norm  is  absolute,  i.e.,  for  every  e  e  ,  s  I — ►  (e(j)s(j))  is  an  isometry.  We 

further  assume  that  the  'unit*  sequences  e^  ,  j  e  J  ,  given  by 

e-'(i)  !-  5^,  all  i,j, 

form  a  basis  for  S  ,  i.e.,  the  truncation  projectr.  p  given  by 


(  y(j>  lien 
(pRy)(j)  «»  < 

l  0  otherwise 

converges  strongly  to  1  as  the  finite  interval  K  approaches  J  •  Then  the  continuous 
dual  S*  of  S  can  be  identified  with  the  sequence  space 

{f  e  RJ  :  Ifl  i*  suPsegf  s/I  si  <  "  }  , 

and  the  norm  on  S*  is  again  absolute.  In  particular, 

I f I  *  I s |  <  lfl*lsl  ,  all  f  6  S*,  s  e  S  . 

Let  A  e  RIXJ  for  some  finite,  infinite  or  biinfinite  integer  intervals  I  and  J 

and  assume  that  A(i,*)  e  S  ,  all  i  .  Then  we  can  identify  A  with  the  linear  map 

*  T 

S  - *  R  s  f  t— ♦  Af  . 

We  are  interested  in  understanding  the  range  of  this  map  under  the  assumption  that  A  is 
SR. 


Theorem  1 .  Let  I ,  J  be  finite,  infinite  or  biinfinite  integer  intervals,  and  let 

S  C  RJ  be  a  normed  linear  space  with  absolute  norm  and  with  ( e^ ) as  a  basis.  If 

Ixj  * 

A  e  R  is  SR,  has  its  rows  in  S  ,  and  carries  some  x  e  S  to  the  strictly 

alternating  sequence  u  :»  Ax  ,  then  the  range  of  A  contains  the  Banach  space 


{v  e  R 


Ivl 


supiel I v( i)/u( i) |  <  “  } 


More  explicitly,  for  every  v  e  Su  there  exists  yv  e  S  so  that  Ayv  =  v  and 

ly  I*  <  Ivl  Ixl*  . 

V  u 


Proof,  we  first  consider  the  case  that  I  is  finite.  Since  S_(u)  =  |l|-1  ,  we  claim 
that  A  has  full  rank  |l|  and  is  therefore  onto.  Indeed,  by  induction,  we  may  assume 

that  A  has  rank  at  least  1 1 1 — 1  .  If  now  rank  A  *  1 1 1  —  1  ,  then  there  would  be,  up  to 

I  .  .  .  # 

scalar  multiples,  a  unique  z  e  R  \101  for  which  z  A  «  0  .  Then  the  sign  regularity  of 

A  would  imply  that  z  must  alternate,  i.e.,  z(i)z(i+1)  <  0  ,  all  i  .  Therefore  0  =  z*Ax 

* 

-  z  u  ,  and  strict  alternation  of  u  would  then  imply  that  z  =  0  ,  a  contradiction. 

It  follows  that  every  v  e  R1  gives  rise  to  a  linear  functional  Fv  defined  on  the 
f inite-dimen8ional  linear  space 


4 


R  s-  span(x(i,‘)) 


tei 


by  the  rule 


£  a  A(i,»)  t— *■  a  v  . 
iel 

In  view  of  the  Hahn-Banach  Theorem,  we  can  therefore  conclude  the  existence  of  yv  e  S 
*  * 

with  Ayv  -  v  and  *  Ivl^lxl  once  we  prove  that 

IF  «  <  Ivl  Ixl*  . 

V  u 

It  Is  sufficient  to  consider  only  finite  J  •  For,  an  infinite  J  can  always  be 
approached  by  finite  intervals  K  ,  and 

A(i,->  -  linW?XMi'*>  '  a11  1  ' 

by  assumption*  Therefore,  for  all  sufficiently  large  intervals  K  ,  the  rule 

£  V’jjAU,*)  *—  a*v 


Fv(r) 


rfiR 


<  sup  inf 
r€R  Ay-v 


reR 

)  y  1  *  I  r  I 

Irl 


and 

lira 

IF*!  - 
K+J  v 

IF  1  . 

V 

of 

Fv 

to  all  of 

S  ,  Then 

« 

sup 

reR 

Ainf  IIH 
Ay-v 

sup 

reR 

'V 

*  If  1 

ly  1* 

'v 

Irl 

with  the  third  inequality  holding  since  the  norm  on  S  is  absolute  and  the  second  since 
Ayy  «  v  .  Now  there  is  no  restriction  in  assuming  that,  for  each  j  e  J  ,  we  can  find 
reR  with  r(j)  ¥  0  (since,  in  the  contrary  case,  the  entire  column  j  of  A  would  be 
zero  and  we  could  simply  omit  the  index  from  consideration).  We  therefore  conclude  from  the 
above  string  of  inequalities  that 


with 


IF  I 

v 


* 

sup  inf  lyl 
wew  Ay-v  W 


W  1-  {r  e  R  !  Irl  <  1,  r(  j)  ?  0  ,  all  j} 

and  1*1  the  dual  norm  to  1*1  ,  i.e., 

w  w 

Ifl*  |f|*|w|  . 
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Now 


while 


.  .*  -  X  8  .  .* 

sup  Ixl  <  sup  J — i — ‘ — L  -  Ixl 
„  w  „„  i  sl 

wew  ses 


inf  lyl 

w 

Ay-v 


IF  I 

v  w 


is  the  norm  of  Fy  with  respect  to  the  I'l^-norm.  It  is  therefore  sufficient  to  prove 
that,  for  any  positive  weight  function  w  , 

inf  lyl*  <  I x I  *  I vl 

w  w  u 

Ay-v 

For  this,  let  yu  be  an  extreme  point  of  the  nonempty,  closed,  bounded,  convex  set 

*  *  * 

C  :=  {yes  :  Ay  -  u  ,  lyl  «  IF  I  }  . 

w  u  w 

We  claim  that  K  :»  supp  yu  contains  exactly  |I|  points.  Indeed,  if  supp  yu  were  to 

contain  a  set  L  with  |L|  «  |l|+1  ,  then  we  could  find  z  e  S  \{0}  with  Az  »  0  and 

*  • 

z^  =  0  .  But  then  ly^+ezl^  -  E  |y^(j)  +  ez<j)lw(j)  "  lyulw  +  e  E  signtyt  j)]  z(  j)w(  j) 
for  all  sufficiently  small  positive  or  negative  c  ,  while  Aly^+ez)  -  u  ,  therefore,  by 
the  minimality  of  lyul  ,  E  sign(y( j ) ] z< j)w< j)  -  0  and  yu  would  not  be  an  extreme  point 
of  C  .  With  this,  (yu)K  is  a  solution  of  the  system  Aj  Ry  =  u  ,  and  our  earlier 
argument  implies  that  Aj  R  is  invertible,  and,  in  particular,  |K|  =  |I|  .  We  therefore 
conclude  that 


inf  lyl  < 
.  w 

Ay-v 


-1 


-  W2WAi,K)"1,k'i,v(i>l  |w(k)l 

<  ,v,u  WieIl^i,Kr'(k^il|“<i>l  |w(k,| 


while 

Wieil^i,K^1tk'i,,,u(i)l  |w(k)| 


EkeK|Eiei(Ai,^’1(k'i)u<i)l  |w(k)l 


ly  I  <  Ixl 
u  w  w 


by  the  sign  regularity  of  A  (which  gives  that  must  be  checkerboard) ,  the 

alternation  of  u  ,  and  the  minimality  of  yu  . 

* 

This  establishes  the  existence  of  y  e  S  with 

v 
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w* 

i 

I 

t 

I 

*  * 

i  Ay  -  v  and  ly  I  <  Ivl  Ixl 

V  V  u 

for  finite  I  •  From  this,  we  obtain  the  result  for  nonfinite  I  by  considering  all  finite 

i  l  * 

1  integer  intervals  L  contained  in  I  .  For  each  such  L  ,  we  can  find  y?  6  S  with 

i  Ay1  -  v.  and  lyLl  <  Iv  I  Ixl  <  Ivl  Ixl  . 

I  V  L  v  L  u  u 

Therefore,  for  some  increasing  sequence  (L)  converging  to  I  ,  the  corresponding  sequence 

L  #  *  *  * 

(y  )  converges  weak  to  some  yes  •  But  then  also  ly  I  <  Ivl  Ixl  and 
Jv  Jv  Jvu 


jxj 

As  a  special  case,  consider  the  SR  matrix  A  e  R  to  carry  to  1^(1)  .  Its 

rows  must  then  be  in  I^(J)  ,  a  sequence  space  with  absolute  norm  and  (e3)  as  a  basis. 

At  the  same  time,  Su  =  1^(1)  provided  u  alternates  uniformly,  i.e.,  u(i)u(i+1)  <  0  , 

all  i  ,  and  inf  |u(i)|  >0  •  We  therefore  have  the  following 


I 


Corollary  1.  _I£  I  and  J  are  finite,  infinite  or  bllnflnlte  integer  Intervals  and 

TwJ 

A  e  R  carries  to  tm(l)  in  such  a  wav  that,  for  some  x  e  ,  u  im  Ax 

uniformly  alternates,  then  A  is  onto. 

Remark.  This  corollary  establishes  the  full  generalization  of  Micchelli's  conjecture. 
The  theorem  even  shows  that  the  solution  y  of  Ay  “  v  may  be  chosen  bounded  in  terms 
of  v  ,  i.e.,  lyl^  <  klvl^  with  k  :«  supi ( j |x( i)/u( j) |  independent  of  v  ,  and  also 
demonstrates  all  this  without  the  assumption  that  A  is  1-1. 


! 

* 

! 

i 

i 


Xx  J 

As  a  second  special  case,  consider  the  SR  matrix  A  €  R  to  have  all  its  rows  in 

■  Cq(J),  another  sequence  space  with  absolute  norm  and  (e^)  as  a  basis e  Then 

*  * 

1*1  «  1*1.  -  1*1  with  e(j)  -  1  ,  all  j  ,  i.e#,  we  are  in  the  special  situation  to 

1  e 

which  we  reduced  the  proof  of  the  theorem#  We  now  know  from  that  proof  that  we  can  choose 
for  each  finite  interval  L  in  I  ,  a  subset  K  of  J  with  |K|  ■  |L|  so  that 
I ( A^  *  min  t 1  :  Ay  -  u  on  L  } 
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■nil  a  irfhilK  mm 
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Then  we  extend  (A.  „)-1  to  CL  e  RJXl  by  talcing  its  values  to  be  zero  off  KxL  .  For 

Li  /  is. 

each  i  e  I  ,  is  in  su  •  The  above  argument  therefore  shows  that 

ICL  e1l1  <  Ix^le1!  -  Ixl | u ( i) |  . 

We  can  therefore  choose  a  sequence  (L)  and  a  corresponding  sequence  (K)  so  that,  for 

each  i  ,  CLe^  converges  weak*  to  some  yA  6  Z^lJ)  .  This  means  that 

•  Li  *  i 

for  all  a  e  c  (J)  ,  lim  a  C  e  =  ay 

U  Ii1*  1 

and  so  in  particular 

(  eX( r)  ,  ieLj 

for  all  r  ,  (Ay^jfr)  »  lim  a(c  e1)(r)  =  lim  i  >  =  eL(r)  . 

L+I  L+I  l  0  ,  i^L  ) 

This  shows  that  the  matrix  C  given  by  C(j,i)  :=  y1(j),  all  (j,i)  e  Jxl  ,  is  the 
pointwise  limit  of  the  sequence  (C^)  .  It  is  a  right  inverse  of  A  and  it  satisfies 

UC( *  ,i)l  <  Ixl^/luti)  |  ,  all  i  . 

This  proves 

Corollary  2.  Let  I  and  J  be  finite,  infinite  or  biinfinite  integer  intervals.  If 
Ixj 

A  e  R  is  SR,  has  its  rows  in  cQ(J)  ,  and  carries  some  x  e  ^(J)  to  the  strictly 

alternating  sequence  u  :=  Ax  ,  then  there  exists  a  sequence  (L)  of  index  intervals 

converging  to  I  and  a  corresponding  sequence  (K)  of  index  sets  so  that  (al  j.)  * 

Jx  X 

exists  and  converges  pointwise  to  a  matrix  C  S  R  which  carries  Su  to  1^(1)  and 
satisfies  AC  =  1  (as  maps,  hence  as  matrices). 

2.  The  inverse  of  a  SR  matrix.  In  this  section,  we  assume  that  the  SR  matrix 

Ixj 

A  e  R  is  also  1-1#  as  a  map  on  ,  in  addition  to  having  a  uniformly  alternating 

sequence  in  its  range*  We  then  know  that  A  is  1-1  and  onto#  hence  invertible,  with  A**^ 

JxI 

again  (representable  as)  a  matrix,  from  R  ,  which  carries  Jl^d)  onto  f„(J)  • 

I  TXT 

Let  now  D  e  R  be  the  diagonal  matrix  whose  diagonal  entries  are  alternately  1 
and  -1  .  Specifically, 

(DIy)(i)  =  (->iy(i)  ,  all  i  €  I,  all  y  e  R1 
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if  I  is  an  interval  (as  we  assume).  It  is  well  known  that,  for  finite  I  and  J  ,  the 


matrix  D^A  'd*  is  again  SR.  In  addition,  if  A  is  TP,  then  DJA  or  its  negative  is 

also  TP.  We  prove  the  first  statement  for  arbitrary  I  and  J  ,  and  prove  the  second 

statement  under  the  additional  assumption  that  the  columns  of  A  are  in  cQ  or  else 
that  I  equals  J  and  is  not  biinfinite,  i.e.,  has  a  first  or  last  entry. 

JXJ 

Proposition  1 •  _I£  A  e  R  maps  i  ( J )  _to  Z^ ( I )  and  is  1-1  and  onto,  and  maps 

Cq(J)  to  Cq(I)  ,  then  A  TP  implies  that  D^A'^D1  or  its  negative  is  TP. 

Proof.  We  know  from  Corollary  2  to  Theorem  1  that,  under  the  given  assumptions, 

A-1  is  the  pointwise  limit  of  certain  matrices  (CL)*  as  the  index  interval  L  converges 

to  I  .  The  matrix  C1,  equals  (A^  j^)  1  =  (A^  ^  )-1  on  Kxi.  and  vanishes  off  KxL  . 

Here  L  is  an  interval,  but  K  is  only  an  index  set,  K  =  {kj,...fk  }  ,  say,  with 

k1  <  ...  <  .  For  such  K  ,  we  define  the  diagonal  matrix  DK  by 

C  C)Ky)  ( kA )  =  (-)1y(ki)  ,  i=1,...,r,  all  yeRK  . 

Then  DL(AKjL)""1DK  is  TP  since  is.  Now  every  i  in  I  must  eventually  be  in  all 

K's  since  in  the  contrary  case  the  i-th  row  of  (CL>*  would  be  zero  for  infinitely  many 

L  ,  hence  A-1(i,«)  =  0  ,  which  is  nonsense.  Thus,  for  any  finite  intervals  M  and  N  , 

(A_1)M,N  is  the  P°intwise  limit  of  (a^)"^  as  t  M  ,  with  ^,ND"( Vi/m.i/  TP 

for  some  e  e  {-1,1}  .  This  implies  that  e„  „  is  independent  of  M  and  N  ,  and  so 
M,  N  M, N 

D**(A  ) m r nDN  or  ^ts  negative  is  TP.  But  since  M  and  N  are  arbitrary  finite  intervals/ 

this  implies  that  D^A^D1  or  its  negative  is  TP.  |  |  ] 

We  believe  the  assumption  that  A  map  Cq(J)  to  Cq(I)  to  be  unnecessary  for  the 

conclusion  that  D**a“^D*  or  its  negative  is  TP.  More  precisely,  we  conjecture  that  a 

Jx  J 

ln- invertible  TP  matrix  A  e  R  must  map  c0(J)  to  cQ(I)  .  Without  this  assumption, 

we  have  no  way  of  approximating  A-1  by  inverses  of  certain  finite  submatrices  of  A  ,  and 
will  have  to  prove  by  some  other  means  that  D^a’^D1  is  SR  in  case  A  is  SR. 
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IXJ 

Theorem  2.  Let  I  ,  J  be  finite,  infinite  or  bllnflnlte  intervals.  If  A  e  R  i* 

SR  and  Invertible  aa  a  map  from  4^1 J)  to  4^(1)  ,  then  DJA  'd1  Is  also  SR. 

Proof.  In  outline,  the  proof  is  as  follows.  By  well  known  results,  it  is  sufficient 
to  prove  that  DJA_,DI  is  variation  diminishing,  i.e., 

s"(dJA~1DIz)  <  S  (z)  ,  all  z  , 
and  this  is  equivalent  to  the  assertion  that 

u  -  Ax  implies  S'(DJx)  <_  S_(DIu)  . 

This,  in  turn,  follows  by  a  smoothing  argument  from  the  assertion  that 

u  »  Ax  and  u,  x  nowhere  zero  implies  S+(D^x)  _£  S+(DIu)  , 

and,  finally,  this  last  statement  follows,  as  we  will  show,  from  the  assertion  that 

u  =  Ax  and  u  nowhere  zero  implies  x  vanishes  at  most  S+(DIu)  times. 

We  begin  the  detailed  argument  with  a  proof  of  this  last  assertion  and  for  this  start 
with  the  following 

1.  if  Be  RIXJ  is  1-1.  then  B_  is  still  1-1  but  not  onto. 

-  ■ 

Proof.  Since  B  is  1-1,  the  sequence  B(»,j)  cannot  be  in  the  range  of  « 

hence  Bj  is  not  onto.  On  the  other  hand,  if  '  t*len'  extending  x 

and  y  to  all  of  J  by  setting  them  equal  to  0  at  j  gives  Bx  »  By  ,  hence 
x  »  y  .  I  I  I 

Corollary,  ^f  u  :=  Ax  uniformly  alternates,  then  x  vanishes  nowhere. 

Proof.  If  x  were  to  vanish  at  j  .  then  the  SR  matrix  would  carry  the 

bounded  sequence  x^  to  the  uniformly  alternating  sequence  u  and  Corollary  1  to  Theorem 
1  would  give  that  is  onto,  while  A  is  1-1  by  assumption,  hence  Ax(jyj  is  not 

onto  by  the  Lemma . | | | 
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Next,  we  strengthen  this  corollary  as  follows 


Proposition  2.  Suppose  u  =  Ax  satisfies  inf^|u(i)|  >  0  and  S+(D1u)  -  k  ,  while 
Xj^  **  0  for  some  L  with  |L|  »  k  .  Let  K  {i  e  1  :  u(i)u(i+1)  >0}  .  Then,  the 
matrix  C  :=  is  again  SR,  1-1  and  onto. 


Proof.  Since  u  never  vanishes  and  S+(DIu)  -  k  ,  therefore  |K|  -  k  and  the 
subsequence  u^  of  u  alternates  uniformly.  In  addition,  Cx^  ■  u^K  and  C  is  SR. 
Therefore  C  is  onto  by  Corollary  1  to  Theorem  1 . 

To  prove  that  C  is  1-1,  let  Cz  «  0  for  some  z  e  t  (\L)  ,  and  extend  z  to 
z  e  i^tJ)  by  z^  =  0  .  Set  y  !“  Az  .  Then  y^  “  0. 

Since  C  is  onto,  we  can  find,  for  each  j  e  l  ,  a  bounded  solution  x3  to  the 
problem 


x3  -  0  ,  Ax3  -  Ae3  off  X  , 
L 


with  e3(i)  {  ,  all  i,j,  as  before.  Set 


F:R 


tm(I)  :  al — *  I  a.  (e3  -  x-*) 

jeL  3 


Then  Fa  =  o  on  L  while  AFa  *  Z  a.(Ae3  -  Ax3 )  vanishes  off  K  .  Therefore, 

3eL  j 

AFa  “  0  on  K  implies  AFa  =  0  and  so,  A  being  1-1,  we  get  Fa  *  0  and,  in 


particular,  a  =  (Fa)r  «  0  .  This  shows  that 

L 


R  .•  a  | — ►  (AFa) 


is  1-1,  hence  onto  since  |L|  «  |K|  . 

It  follows  that  we  can  choose  a  so  that  AFa  -  y  on  K  .  But  then  z'  z  -  Fa 

satisfies 


f  y  -  y  on  x| 

Az'  “  y  -  AFa  -  \  S  -  0 

\  0  -  0  off  K  J 

and  so,  A  being  1-1,  we  have  z'  “  0  ,  therefore  0  “  z'  ■  0  -  a  ,  i.e  ,  a  •  0  and  so, 

Xi 

finally,  z  -  z^  -  (z  -  Fa)NL  *  z^  =0  -  III 


! 


Remark.  The  argument  just  given  shows 
B  is  1-1  and  can  be  partitioned  as 

B 

in  such  a  way  that  is  onto  while  Bj2 

also  1-1. 

Corollary.  .If  u  =  Ax  e  1^(1)  with  inf^Juli))  >  0  and  S+(D1u)  =  k  ,  then  x  has 
at  most  k  zero  entries. 

Proof.  Let  Xj^  -  0  for  some  L  with  |L|  =  k  .  Setting  again 
K  :=  {i  e  I  s  u(i)u(i+1)  >0}  ,  we  know  from  Proposition  2  that  C  a\jc,\l  is 

while  it  obviously  carries  x^L  to  the  uniformly  alternating  sequence  u^  and  is  SR. 
Therefore,  by  the  Corollary  to  Lemma  1,  x  does  not  vanish  off  L  .III 

.  2.  .If  u  =  Ax  with  infjjua if  >  0  and  S+<DIu)  *  k  ,  then  S+(DJx>  <  k  . 

Proof.  We  first  show  that  we  may  assume  that  x  vanishes  nowhere.  For,  if  this  is  not 

the  case,  then  we  replace  each  zero  entry  of  x  by  e  or  -e  in  such  a  way  that  the 

resulting  sequence  xE  satisfies  S+(DJxe)  =  S+(DJx)  .  This  changes  u  =  Ax  to 
u6  :=  Ax6  =  u  +  v  with  lvlm  <  IAI|e|  .  But  since  infjjuf  i)  I  >  0  ,  we  can  choose  e  >  0 
so  small  that  again  inf^|u£(i)|  >  0  ,  while  S  (DIuE)  =  S  (D  u) 

Next  we  produce  a  SR  1-1  onto  matrix  C  and  a  sequence  z  with  as  many  zeros  as 

DJx  has  sign  changes  and  with  Cz  =  u  .  For  this,  consider  the  matrix  E^a)  which 

differs  from  the  identity  only  in  that  it  has  an  a  in  position  (1+1,1)  .  This  matrix  is 
TP  for  nonnegative  a  and  carries  the  sequence  z  to  itself  except  that  the  (i+1)st  entry 

is  changed  to  z(i+1)  +  az(i)  .  Consequently,  Ei(a)  is  invertible,  with  Ei(-a)  its 

inverse. 


the  following  general  fact:  If  the  linear  map 


is  square  of  finite  order  k  ,  then  B11  is 
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Now  let  :=  x(i+1)/x(i)  .  If  x(i)x(i+1)  >  0  ,  then  >  0  and  y  E^(-r^)x 

equals  x  except  for  a  zero  in  entry  i+1  .  Hence,  if  i1  <  ...  <  in  are  all  in 
K  j“  {i  e  I  s  x(i)x(i+1)  >  0}  ,  then  the  matrix 


E^  (  — r,  )  ..»  E,  t-r, 
11  11  xn  *n 


carries  x  to  a  sequence  which  vanishes  at  i^  +  1,  ...»  in+1  .  while 

-1 


B 


E.  (r.  )  ...  E.  (r.  ) 
i  i  i.  i. 

n  n  11 


is  TP,  1-1  and  onto,  hence  AB-1  is  again  SR,  1-1  and  onto.  Since  AB-1(Bx)  =*  u  ,  we  now 


conclude  from  the  Corollary  to  Proposition  2  that  n  <  k  .  This  proves  the  lemma  in  view  of 
the  fact  that  s+(DJx)  -  |K|  ,  since  x  vanishes  nowhere.  I  I  I 


3.  If  Ax  «  u  ,  then  S'fD^x)  <  S'tD^u). 


Proof.  There  is  nothing  to  prove  unless  s“(DIu)  <  «>  .  In  that  case,  we  choose 
sign[u(i)]  e  {-1,1}  in  such  a  way  that  S+( D1 ( sign! u( • ) ] ) )  *  S_(DIu)  and  then  set 


u£(i) 


|  e  sign(u(i)J  if  lu(i)|  <  c 

l  u(i)  otherwise 


Then  S+(DIu£)  =  S  (D^u)  and  so,  using  the  boundedness  of  A-'  , 

S  (DJX)  <  lim  S+(DJ(A  'u6))  <  S+(DIuE)  »  S  (DIu)  , 


by  Lemma  2.  |  |  | 


With  this,  the  proof  of  Theorem  2  is  apparent.  For  we  now  conclude  from  Lemma  3  that 
s”(DJA_1DIz)  <  s”(z)  ,  all  z  ,  and  therefore  every  finite  submatrix  of  DJA-1DI  is 
variation-diminishing.  Hence,  by  Karlin  (3ip.222J,  D^A^D1  is  SR.  ||| 

Corollary.  lt_  I  -  J  i3  only  Infinite  (and  not  bllnfinitel,  then  A  TP  implies 
DA-1D  TP. 

Proof.  Assume  without  loss  that  I  •  {1,2,3,...}  and  consider  the  matrix  B  e  rKxK 


1  3 
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